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Abstract The quark combination mechanism of hadron production is applied to nucleus-nucleus collisions at 
the CERN Super Proton Synchrotron (SPS) and BNL Alternating Gradient Synchrotron (AGS). The rapidity 
spectra of identified hadrons and their spectrum widths are studied. The data of ty~ , K ± , <f}, A, A, H - , and H + 
at 80 and 40 AGeV, in particular at 30 and 20 AGeV where the onset of deconfinement is suggested to happen, 
are consistently described by the quark combination model. However at AGS 11.6 AGeV below the onset 
the spectra of ir , K and A can not be simultaneously explained, indicating the disappearance of intrinsic 
correlation of their production in the constituent quark level. The collision-energy dependence of the rapidity 
spectrum widths of constituent quarks and the strangeness of the hot and dense quark matter produced in 
heavy ion collisions are obtained and discussed. 
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1 Introduction 



The production of quark gluon plasma (QGP) and 
its properties are hot topics in relativistic heavy ion 
collisions. A huge number of possible QGP signals 
were proposed and measured, and many unexpected 
novel phenomena were observed at RHIC and SPS 
[lJ-Q. These experimental data greatly contribute to 
the identification of QGP and the understanding of its 
properties and hadronization from different aspects. 
Especially, there are a class of phenomena that are 
of particular interest, i.e. the abnormally high ratio 
of baryons to mesons and the quark number scaling 
of hadron elliptic flows in the intermediate p T range, 
etc@, 0]. They reveal the novel features of hadron 
production in relativistic heavy ion collisions. 

In quark combination/coalescence scenario, 
hadrons are combined from quarks and antiquarks, 
i.e., a quark-antiquark pair merges into a meson and 
three quarks into a baryon. The production difference 
between baryon and meson mainly results from their 
different constituent quark numbers. It is shown that 
such a simple quark number counting can naturally 
explain those striking features of hadron production 
observed at RHIC[8l4l0j . while the fragmentation 



mechanism can not. 

The highlights at RHIC are mainly of hadron pro- 
duction in transverse direction where the quark com- 
bination scenario mostly flashes. In fact, the longitu- 
dinal rapidity distribution of hadrons is also a good 
tool for tes ting the hadronization mechanism. In pre- 
vious work 11 , li| , we have used the quark combina- 
tion model to successfully describe the rapidity spec- 
tra of various hadrons at RHIC y/s NN = 200 GeV 
and top SPS E beam = 158 AGeV. At other collision 
energies where the QGP may be produced, e.g., at 
lower SPS and higher LHC energies, does the quark 
combination mechanism still work well? The Beam 
Energy Scan at RHIC and the NA49 Collaboration 
have provided abundant data on hadron production 
in the energy region from 20 GeV to 6 GeV. In this 
paper, we extend the quark combination model to 
systematically study the rapidity distributions of var- 
ious identified hadrons in heavy ion collisions at SPS 
E beam = 80, 40, 30, 20 AGeV and AGS E beam = 11.6 
AGeV and test the applicability of the quark combi- 
nation mechanism. 
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2 A brief introduction to the quark 
combination model 

The quark combination model deals with how 
quarks and antiquarks covert to color-siglet hadrons 
as the partonic matter evolves to the interface of 
hadronization. The basic idea is to put all the quarks 
and antiquarks line up in a one-dimensional order in 
phase space, e.g., in rapidity, and let them combine 
into primary hadrons one by one following a combi- 
nation rule based on the QCD and near-correlation in 
phase space requirements. See Sec. II of Ref.[12l] for 
detailed description of such a rule. Here, we consider 
only the production of SU(3) ground states, i.e. 36- 
plets of mesons and 56-plets of baryons. The flavor 
SU(3) symmetry with strangeness suppression in the 
yields of initially produced hadrons is fulfilled in the 
model. The decay of short-life hadrons is systemati- 
cally taken into account to get the spectra comparable 
to the data. The model has reproduced the experi- 
mental data for hadron multiplicity ratios, momen- 
tum distributions and the elliptic flows of identified 
hadrons, etc., in heavy ion collisions at RHIC and 



issue 



top SPS energies [P2]4l6j], an d addressed the entropy 



17| and the exotic state production[l 



3 Rapidity spectra of constituent 
quarks 

The rapidity spectra of constituent quarks just 
before hadronization are needed as the input of the 
quark combination model. Considering that the colli- 
sion energies studied here are much lower than RHIC 
energies, and in particular 30-20 AGeV is the possible 
region for the onset of deconfinement, it is not sure 
whether the hot and dense quark matter is exactly 
produced at these energies, so applying a model or 
theory for the evolution of the hot and dense quark 
matter, e.g., relativistic hydrodynamics, to get the 
quark spectra before hadronization may be uneco- 
nomic or infeasible. In this paper, assuming the 
hot quark matter has been created, we parameterize 
the rapidity distribution of constituent quarks and 
extract the parameter values from the experimental 



data of final state hadrons. The rapidity distribu- 
tion of newborn quarks is taken to be a Gaussian-like 
form, i.e., 



dN q 
dy 



NJ q {y) 



A 



c 



(i) 



Here C = exp[—\y bearn \ a /2a 2 } which means the con- 
stituent quarks that form hadrons via combina- 
tion are within [— i/t eam ,^ elim ] in the center-of-mass 
frame. A is the normalization factor, satisfying A = 



J yb — (e 



-C)dy. a and a are the parameters 



depicting the shape of the spectrum. N q denotes the 
total number of newborn quarks with type q in the 
full rapidity region. For the net quarks coming from 
the colliding nuclei, their evolution is generally dif- 
ferent from the newborn quarks due to complicated 
collision transparency 19(. We fix the rapidity spec- 
trum of net-quarks before hadronization by the data 
of the rapidity distribution of fnet-'lprotonsfiol]. and 
the results are shown in Fig. [TJ 
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Fig. 1. The rapidity distributions of net-quarks. 

4 Hadronic rapidity spectra at SPS 
and AGS 

In this section, we first use the quark combination 
model to systematically study hadron rapidity spec- 
tra and their widths at SPS energies. Then we present 
the energy dependence of the strangeness and spec- 
trum width of constituent quarks. Finally, we show 
the results at AGS 11. 6 AGeV. 

4.1 Rapidity spectra of hadrons at SPS ener- 
gies 
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Fig. 2. Rapidity distributions of identified hadrons in central Pb+Pb collisions at Eteam = 80, 40, 30, 20 
AGeV. The symbols are the experimental data [2l| - [25| and lines are the calculated results. The open symbols 
are reflection of measured data (solid symbols) against midrapidity. The results of <j> & t 30, 20 AGeV are 
scaled by proper constant factors for comparing their shapes with the data. 



The calculated rapidity distributions of identified 
hadrons in central Pb+Pb collisions at E beam = 80, 
40, 30, 20 AGeV are shown in Fig. [2j The values 
of parameters for newborn quarks are listed in Table 
[U and x 2 /ndf is presented as the quality of global 
fitting. Note that the experimental data of pions be- 
yond Ubeam n ot included in x 2 / n df calculation, 
because the limiting fragmentation behavior around 



Vbeam. becomes prominent which is beyond the study 
of this paper by quark combination. The results show 
the quark combination model can reproduce the ex- 
perimental data of various identified hadrons except 
<f> at 30, 20 AGeV. The yields of at 30, 20 AGeV 
deviate from the data, but their sharps are still in 
agreement with the data after scaling by proper con- 
stant factors. 



Table 1. Parameters of newborn light and strange quarks and % 2 /ndf '. 



energy 


a u/d 


a s 




c s 


N u/d 


N s 


X 2 /ndf 


80AGeV 


1.90 


2.35 


1.19 


1.40 


275.2 


159.6 


0.98 


40AGeV 


1.90 


2.25 


1.12 


1.40 


171.3 


116.5 


1.01 


30AGeV 


1.80 


2.05 


1.20 


1.16 


139.0 


101.5 


0.96 


20AGeV 


1.80 


1.95 


1.10 


0.85 


105.5 


77.0 


1.90 
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The deviation of 4> yields at 30 and 20 AGe V 
is related to the pronounced rescattering effect. It 
is shown that at higher SPS and RHIC energies 
the production of <p mainly comes from the con- 
tribution of partonic phase, i.e., the directly pro- 
duced <j) by hadronization while at lower SPS energies 
kaon coalescence may be the dominated production 
mechanism 25]. In addition, the directly produced <j) 
after hadronization will possibly suffer the destroy by 
the scattering of the daughter kaons with other pro- 
duced hadrons. The absence of these two effects at 
hadronic stage in our calculations is the main reason 
for the deviation in <f> yields. 

4.2 Widths of rapidity spectra for hadrons 

From the above experimental data one can see the 
widths of rapidity spectra for different hadron species 
are generally different. This difference, in other word 
the correlation of longitudinal hadron production, can 
be used to quantitatively test various hadronization 
models. In quark combination mechanism, the rapid- 
ity distribution of a specific hadron is the convolu- 
tion of the rapidity spectra of its constituent quarks 
and combination probability function (denoted by the 
combination rule in our model). Since the rapidity 
distributions of different-flavor quarks are different 
(see Tableland Fig. [T]), in particular the difference 
between newborn quarks and net-quarks, the shapes 
of rapidity spectra of hadrons with different quark 
flavor components are generally different and corre- 
lated with each other by constituent quarks. Here, 
we calculate the spectrum widths of various hadrons 
to clarify this feature. 

Considering that the rapidity region covered by 
the data for different hadron species or at different 
collision energies are not the same, we define the vari- 



ance of rapidity distribution for a specified hadron in 
a finite rapidity region limited by the discrete exper- 
imental data, 



E Vi dy 

(y 2 ) (L) = — n ■ 



(2) 



•hj 



Here n is the number of experimental data; yi and 
dNi/dy are the rapidity position and the correspond- 
ing yield density measured experimentally, respec- 
tively. Replacing dNi/dy with the model results, we 
can give the (y 2 }^ \ mo dei and compare it with the ex- 
perimental value (y 2 ) (L) |data to test the applicability 
of the model without any arbitrariness caused by the 
rapidity region where the experimental data have not 
covered yet. We further extrapolate {y 2 ) {L) \data to the 
full rapidity region [~y beam ,ybea m ] via the relation 



(y 2 ) (F) \ 



(y 



2\(L)\ 



(y 



2\(i) 



-(y 



2\(F)| 



(3) 



where (y 2 )^ \ mo dei is the variance calculated by the 
model in the full rapidity region. The degree of agree- 



ment between (y 



2\(F)| 



and experimental data 



(y 



2\(F)| 



i ata still represents the original description 



ability of the model, and the < y 2 > of different 
hadron species can be directly compared and their en- 
ergy dependence is recovered. The spectrum widths 
of various hadrons, i.e. D(yY F) = sj (y 2 ) (F) , are cal- 
culated and the results are shown in Fig|3] (158AGeV 
is also included). One can see that the D(y) of var- 
ious hadrons given by the model are obviously dis- 
tinguished. The agreement between the data and 
the calculated results is the support of the theo- 
retic (quark recombination mechanism) explanation 
for the widths of hadronic rapidity distributions. 




Fig. 3. Widths of rapidity distributions D(y) (F) at E beam (y beam ) = 158, 80, 40, 30, 20 AGeV (2.91, 2.57, 
2.22, 2.08, 1.88). The right panel includes the results of K + K~ coalescence for <f> production. Model results 
are connected by the solid lines to guide the eye. 
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As stated above, at lower SPS energies, kaon me- 
son coalescence may be an important mechanism for 
(f> production. The right panel of Fig. [3] shows the 

as well as the re- 
(F) \coai- Similar to 



(F)\ 



4's D(y)^\ data and D(y) 
suits for K + K~ coalescence D(ij) 



Ref .[25], considering the ideal case of coalescence of 
two kaons with the same rapidity, we use the mea- 
sured kaon rapidity distributions f K ±{y) to obtain 
the spectrum of <f> by f^ al {y) oc f K +{y)f K -{y) and 
then give the D(y) (F) \ coa i- One can see that at col- 
lision energies above 20 AGeV, D(yY F) \ coal is much 
smaller than the data, which is similar to the results 



in Ref. 25], while our results nearly agree with the 
data. This clearly shows the </> production at these 
energies is dominated by the hadronization. At 20 



AGeV, D(y) 



is nearly equal to the data and 



the model result is also in agreement with the data. It 
suggests (f> production at lowest SPS energy can have 
different explanations, in other words, even though 
kaon coalescence is significant, the <fi directly pro- 
duced from hadronization may be unnegligible. 



4.3 The strangeness and spectrum width of 
constituent quarks 

Let us turn to the extracted rapidity distributions 
of newborn light and strange quarks. Their properties 
can be characterized by two quantities, i.e the ratio 
of strange quark number N a to light quark number 
N u / d called strangeness suppression factor A s and the 
width of rapidity spectrum. 

The left panel of fig. H] shows the strangeness 
suppression factor X s at different collision energies. 
Note that A s defined here is in terms of quark num- 
bers in the full phase space, so the values are sightly 
different from those in terms of midrapidity quark 
number densities in previous Ref. 12|. As the 
comparison with SPS, we also present X s at RHIC 
^/Sjvjv = 200, 62.4 GeV calculated by the model and 
at AGS 11.6 AGeV calculated by counting the num- 
bers of light and strange valance quarks hidden in the 
observed pions, kaons and A which are mostly abun- 
dant hadron species carrying light and strange gradi- 
ents of the system. One can see that the value of A s 
exhibits a peak behavior at lower SPS energies. This 
behavior has been reported by NA49 Collaboration 
as the signal of onset of deconfinement. 
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Fig. 4. The strangeness suppression factor A s and D(y) of newborn light and strange quarks at different 
collision energies. The results are connected by the solid lines to guide the eye. 



The width of rapidity spectrum is characterized 
by D{y) = Jtf 



(y) 2 = y (y 2 )- The right panel of 
Fig. U] shows the D{y) of newborn light and strange 
quarks at different collision energies. The results at 
RHIC y/IJ^ = 200, 62.4 GeV and SPS E beam = 158 
AGeV are included. One can see that with the in- 
creasing of collision energy D{y) of newborn light 
and strange quarks both increase regularly, indicat- 
ing the stronger collective flow formed at higher en- 
ergies. As collision energy is equal to or greater than 
40 AGeV, D(y) of newborn light quarks are always 
smaller than that of strange quarks while at 30 and 



20 AGeV the situation flips. The wider spectrum 
of strange quarks, relative to light quarks, has been 
verified at RHIC both in the longitudinal and trans- 
verse directions 12 , [3, 26 1 , and the explanation is 



that strange quarks acquire stronger collective flow 
during evolution in partonic phase. As the collisions 
energy reduces to 30 and 20 AGeV, the widths of ra- 
pidity distributions of light and strange quarks, see 
Table [1] and Fig. |4l are all quite narrow, which 
means the collective flow formed in partonic phase 
is much smaller than those at higher SPS and RHIC 
energies. Therefore, the partonic bulk matter, even 
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if produced as the indication of our results via still 
active constituent quark degrees of freedom, should 
be in the vicinity of phase boundary, and the ex- 
tracted momentum distributions of quarks keep the 
memory of their original excitation. If thermal equi- 
librium is reached in heavy ion collisions, the quark 
occupation function in momentum space follows as 
exp[— E/T] — exp[— m T cosh(y)/T] in the case of no 
collective flow. Taking the hadronization tempera- 
ture T = 165 MeV and constituent mass m q = 340 
MeV for light quarks and m s — 500 MeV for strange 
quarks, the quark rapidity spectrum is Gaussian form 
and width of light quarks is <r q = 0.6 and strange 
quarks a s — 0.52 due to the heavier mass. The 
tighter spread of strange quarks in rapidity space can 
be qualitatively understood in quark production with 
thermal-like excitation. 

4.4 Results at AGS 11.6 AGeV 

What happens at lower AGS energies? We fur- 
ther use the model to calculate the rapidity distribu- 



tions of various hadrons at 11.6 AGeV. The results 
are shown in Fig. [5] and are compared with the ex- 
perimental data. The values of parameters (a,a,N q ) 
for quark spectra are taken to be (2.1, 0.88,71) for 
newborn light quarks and (2.0, 0.83, 42) for strange 
quarks, respectively. The rapidity distribution of net- 
quarks is extracted from the proton data of E802 Col- 
laboration QUof, and the data of E877 Collabora- 
tion [i^l at 10.8 AGeV are used as the guide of extrap- 
olation of net-quark spectrum in the forward rapidity 
region. One can see that the results for pions and 
kaons are in agreement with the data but the result 
of A can not reproduce the experimental data — the 
width of spectrum given by the model is much wider 
than that of data. This suggests that there is no in- 
trinsic correlation at constituent quark level between 
production of kaons and A at AGS 11.6 AGeV. In ad- 
dition, the rapidity distributions of 0, K°, p, A and 
S~(2 + ) are predicted to be tested by future data. 




5 Summary 

In this paper, we have investigated, with the quark 
combination model, the rapidity distributions of iden- 
tified hadrons and their widths in central A+A colli- 
sions at SPS and AGS energies. Assuming in advance 
the existence of constituent quark degrees of freedom, 
we parameterize the rapidity spectra of quarks be- 
fore hadronization, then test whether such a set of 
light and strange quark spectra can self-consistently 



explain the data of ir~, K ± , (j>, A(A), S~(S ), etc., 
at these energies. The results of hadronic rapidity 
spectra are in agreement with the data at 80 and 40 
AGeV. At 30 and 20 AGeV where the onset of decon- 
fincment is suggested to happen, the model can still 
basically describe the production of various hadrons. 
The study of rapidity-spectrum widths for hadrons, 
particularly for phi via the comparison to the results 
from kaon coalescence in the stage of hadronic rescat- 
terings, clearly show the hadron production at the 
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collision energies above 20 AGeV is dominated by 
the hadronization. The energy dependence of the 
rapidity-spectrum widths of constituent quarks and 
the strangeness of the hot and dense quark mat- 
ter are obtained. It is shown that the strangeness 
peaks around 30 AGeV and below (above) the energy, 
the width of strange quarks becomes to be smaller 
(greater) than that of light quarks. As the collision 
energy decreases to AGS 11.6 AGeV, it is found that 
the production of n ± , K ± and A can not be consis- 
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